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ABSTRACT 

We use a gauge-invariant regularization procedure, called split dimensional regularization, to 
evaluate the quark self-energy E(p) and quark-quark-gluon vertex function A /1 (p',p) in the 
Coulomb gauge, v ' A a = 0. The technique of split dimensional regularization was designed 
to regulate Coulomb-gauge Feynman integrals in non-Abelian theories. The technique which 
is based on two complex regulating parameters, uo and a, is shown to generate a well-defined 
set of Coulomb-gauge integrals. A major component of this project deals with the evaluation 
of four-propagator and five-propagator Coulomb integrals, some of which are nonlocal. It is 
further argued that the standard one-loop BRST identity relating E and A M , should by rights 
be replaced by a more general BRST identity which contains two additional contributions 
from ghost vertex diagrams. Despite the appearance of nonlocal Coulomb integrals, both E 
and A M are local functions which satisfy the appropriate BRST identity. Application of split 
dimensional regularization to two-loop energy integrals is briefly discussed. 
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1 Introduction 



In our first test of the technique of split dimensional regularization, we evaluated the Yang- 
Mills self-energy to one loop and verified the appropriate BRST identity 0. In this article, 
we present a second test of split dimensional regularization, also in the Coulomb gauge, by 
calculating the quark self-energy and quark-quark-gluon vertex functions. Our rationale for 
concentrating on one-loop calculations at this stage is very simple: we want to make sure that 
split dimensional regularization is capable of regulating unambiguously all one-loop integrals 
in the Coulomb gauge, before tackling the infamous divergences at two loops and beyond. 

The Coulomb gauge has been highly successful in both electrostatics and quantum electro- 
dynamics, especially in the treatment of bound-state problems ||. By contrast, its application 
to non-Abelian models has been marred by the persistent appearance of ambiguous Coulomb- 
gauge Feynman integrals. It was Schwinger || who first noticed that the transition from 
a classical Hamiltonian to a quantum Hamiltonian led to instantaneous Coulomb-interaction 
terms. These terms, subsequently labelled Vi + "VVterms by Christ and Lee [|J , are caused by 
the operator ordering problem, and arise whenever we convert a classical system to a quantum 
system. The principal problem with these new Vi + V2-terms is their nonlocality, foreshadowing 
potential difficulties in evaluating Feynman integrals. 

A review of the historical development of the Coulomb gauge in the framework of Yang- 
Mills theory can be found in Ref. and will not be repeated here. Suffice it to say, however, 
that the nonlocal Vi + "VVterms of Christ and Lee lead in general to momentum integrals 
containing two types of divergences: (a) ordinary UV divergences, related to the structure of 
space-time, and (b) new divergences, arising from the integration over the energy variable dqo, 
where d 4 q = dqo d 3 q. Hence the name energy-integrals || ||. It is the divergences from the 
energy-integrals that cause the notorious ambiguities in the Coulomb gauge, especially at two 
and three loops. To compound the issue, standard dimensional regularization is not powerful 
enough to consistently regulate the two types of infinities mentioned above, as first pointed out 
by Doust and Taylor ]7|, ||. 

In 1996, we introduced a radically different method, called split dimensional regularization, 
for regulating Feynman integrals in the Coulomb gauge and applied it to the one-loop gluon 
self-energy IP™ [[[], |9], [10| . The new technique employs two complex-dimensional parameters uj 
and a, replacing the measure d 4 q = dqo d 3 q by 

where the limits o — > (l/2) + , and uo — > (3/2) + are to be taken after all integrations have been 
completed. Application of this procedure to the gluon self-energy Ilf^(p) leads to the following 
results QTJ] : 

(i) IE™ is local and nontransverse in the Coulomb gauge. 

(ii) Ghosts play a significant role, despite the physical nature of the Coulomb gauge. 

(iii) The BRST identity contains a nonvanishing ghost diagram, but is rigorously satisfied. 

Encouraging as this initial result may be, it would be foolhardy to assume that split dimen- 
sional regularization can regulate not only the remaining one-loop integrals in the Coulomb 
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gauge, but two- and three-loop integrals as well. After all, it is well known how annoyingly 
unpredictable noncovariant gauges may be at times! For this reason, we have decided to pro- 
ceed with the investigation of one-loop diagrams, computing both the quark self-energy and 
quark-quark-gluon vertex function. 

The latter is particularly challenging, since it contains for the first time, Coulomb-gauge inte- 
grals with five propagators. 

The paper is organized as follows. In Section 2, we review the Feynman rules in the 
Coulomb gauge, and evaluate the quark self-energy S(p) using split dimensional regulariza- 
tion. The quark-quark-gluon vertex function A M (p',p) is calculated in Section 3. In Section 4, 
we illustrate the technique of split dimensional regularization for more complicated Coulomb- 
gauge integrals: for a massless four-propagator integral containing the instantaneous Coulomb- 
interaction propagators (q 2 )^ 1 and [(q + q) 2 ]^ 1 , and for a massive five-propagator integral with 
propagators (q 2 )" 1 and [(q — k) 2 ] -1 . In Section 5, we examine the ghost contributions to the 
BRST identity and verify the latter. Implications of the technique of split dimensional regular- 
ization for two-loop energy integrals are discussed in Section 6. The highlights of our calculation 
are summarized in Section 7. The Appendix contains a partial list of integrals needed for the 
determination of S(p) and A M (p',p). 

2 Feynman rules and computation of quark self-energy 

2.1 Review of Feynman rules [11,12,13] 

In Yang-Mills theory, with the Lagrangian density 

V = L - —(F* h k h ^) 2 , a = gauge parameter, a -> , (2) 
2a ^ 

where 

Ff = (^-^ % )5 a \ ^ = 0,1,2,3, 
Ff A" b = V • A a , n M = (no, 3) = (l,0),n 2 = 1 , 

L = "^(F^) 2 + (J< + c a F; c )D cb ^ b - igf abc K a cV , 

F% = d,K - d u A; + gf abc A b A^ , 
D ab = 5 ab «9 M + gf abc A^ , 

the noncovariant Coulomb gauge is given by 

V-A a (x) = 0; (3) 

A a is a massless gauge field, g the gauge coupling constant, f abc are group structure constants, 
and a = 1, ...,N 2 — 1, for SU(N); uj a ,u a represent ghost, anti-ghost fields respectively, while K a 
and J a denote external sources. The quantities J a , u & and a) a are anti-commuting. 
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In the Coulomb gauge, the Lagrangian density L' leads to the following gauge boson prop- 
agator PJ: 



G£(?) 



-id 



ab 



(27r) 4 (q 2 + ie) 
while the three-gluon vertex is given by 



- q ■ n(q / ,n i , + q^) 

-q 2 



e > , (4) 



VS(p,?,r) = gf abc (27r) 4 5 4 (p + q + r)[g^(p-q), 



fabc/Q„\4 r4 



+ g*p(q - r)„ + g„ M (r - p)„ 



(5) 



and the scalar ghost propagator by 



; £ab 
pab _ 10 
ghost /o_N4^2 



(2vr) 4 q 

2.2 Quark self-energy 

The unintegrated expression for the quark self-energy in Minkowski space reads (cf. Fig. 1) 



(6) 



^Coul 



(P) 



4i g 2 



d 4 q 



3 J (27r) 4 (q 2 + ie) /j-m + ie 

^ ( n 2 q M q^ — q • n(q M n 1 ' + q"n M 



—7m 

, n„=(l,0,0,0),e>0, (7) 



where m is the quark mass; the generators T^g have already been multiplied out. Expression 
(7) leads to trivial covariant integrals, and to noncovariant Coulomb-gauge integrals involving 
the factor 1/q 2 . By performing a Wick rotation and applying split dimensional regularization 
to all Coulomb integrals, we may derive the following formulas in Euclidean space: 



div 
div 
div 
div 



dq 


Kp- 


q) 2 + m^q 2 




dq ft 


Kp- 


q) 2 + m^q 2 




dq 


q 2 [(p 


— q) 2 + m 2 ]cp 




dq ft^j 



21* , dq 



( 27r )2(a+^) 



- Pi r, i = 1,2,3, 



— i* 



where 



I* = div 



q 2 [(p — q) 2 + m^q 2 3 1J 
d 2uj q f d 2a 



(8a) 
(8b) 
(8c) 
(8d) 



(2tt) 2 ^ 
divergent part of 
T{2-u-a 



(27r) 2CT q 2 (q + p) 2 ' 

T{2 - u - a)(p 2 )^- 2 



(4?r) ; 



(47r) w+CT 
uj -(3/2)+, a ^(1/2)" 



(9) 
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Executing the required integrations in Eq. (7), we obtain 

s Coul (p) = 1 ( ^ _ 4m)r(2 _ u _ ff)i (10) 

with a s = g 2 /47r , g = g/i 2_w_<T , and \x being the mass scale. We see that the result in Eq. (10) 
is strictly local, despite the appearance of nonlocal integrals such as Eq. (8c). 



3 Quark-quark-gluon vertex 

The quark-quark-gluon vertex function A M (p', p) consists of the QED-like quark-quark-gluon 
vertex, r^(p',p), and the non-Abelian quark-quark-gluon vertex, r 2 (p',p): 

A,(p',p) = ri(p',p) + rJ(p', P ) . (ii) 

r* (Fig. 2) is the easier one to compute, containing as it does only the single noncovariant 
factor 1/q 2 . Thus, 

r i/ / \ ig 2 f d 4 q 1p A + m hv(£>- k + m )7a 

^ P ' PJ 6(27r) 4 i q 2 [(p'-q) 2 -m 2 ][(p-q) 2 -m 2 ] 

'nVn CT - q • n(q a -a p + q p n a y ' 



(12) 



where the ie-terms in the covariant propagators have been omitted for clarity. The divergent 
component of r* is, again, amazingly simple: 



r i(P> P) = |^(2 - u - a), u - (3/2)+, a -> (1/2)+ . (13) 

In deriving this answer, we have made use of the formulas listed in the Appendix. Notice in 
particular the nonlocal integrals there. 

Computation of the non-Abelian vertex function r 2 (p',p) in Fig. 3, 

r >' P) = 2W^ q 2 (k-q) 2 [(p + q) 2 -m 2 ] - {g + ^ ^ ~~ ^ n(q p n + q^)]} 
• |g A<J + -=^^[n 2 (k - q) A (k - q)' - n • (k - q)[(k - g ) V + (k - q) CT n A ]]| 
•[-(k + q) A g Mp + (2q-k) M g pA + (2k-g) p g MA ] , k = p' - p , (14) 

is complicated by the presence of the two noncovariant factors 1/q 2 and l/(k — q) 2 , and the 
necessity of having to evaluate for the first time /we-propagator integrals (see Section 4 for 
details). In fact, determination of the pole parts of these five-propagator integrals in the 
Coulomb gauge took as long as the rest of the calculation. 

Despite the presence of explicitly nonlocal integrals in the evaluation of T 2 , Eq. (14), the 
final expression turns out to be strictly local: 

r>',p) = ^7^(2 - a; - (7) , (15) 
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so that 



A„(p',p) 



-Q! s 



3tt 



7 M r(2 - w - a) 



a->(l/2)+ 
cj-*(3/2)+ 



(16) 



4 Procedure for noncovariant four-propagator and five- 
propagator integrals 

In this section, we shall demonstrate our procedure for Coulomb-gauge integrals containing 
precisely two noncovariant factors of the form 1/q 2 and l/(k — q) 2 . We shall first tackle 
the four-propagator integral, then continue with the five-propagator case, consisting of three 
covariant and two noncovariant propagators. 

4.1 The four-propagator integral 

Application of split dimensional regularization to the Euclidean-space integral I, 



yields [|4 



I 



(27r)2(«-Hr) q 2 (q + p)V(q + p) 2 



(2tt) 2w 7 (2tt) 2ct 
■exp{-A[q 2 + 2yq ■ p + yp 2 + xq 2 + 2Gq 4 p 4 + Gp 2 ]} 



;i7) 



with 



— — - - / dx / dG / dy / dAA 3 -"- CT x- c 



(2tt) 2 (' 

•exp{-A[(y-y 2 )p 2 + G(l-G/x)p 2 ]} 



7T 



uj-\-a 



(19) 



H(x, G, y) = y(l - y)p 2 + G(l - G/x)p 2 . (20) 

In order to be able to integrate over y, G and x in Eq. (19), we need to pinpoint all possible 
singularities in parameter space. To this effect, we observe that the integral I diverges for H = 
0, i.e. for the following two cases pj: 



Case 1 : H -> 0, if y = and G = ; 
Case 2 : H — > 0, if y = 1 and G = x . 



(21a) 
(21b) 
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Since H approaches zero linearly in Eq. (20), we may drop the terms proportional to y 2 and 
G 2 in Case 1, so that 

H -> H? } = yp 2 + Gp 2 . (22a) 

Similarly, setting 1 — y = Y, and x — G = t in Case 2, we may drop the terms proportional to 
Y 2 and t 2 , in which case 

H^H( 2) = (l-y)p 2 + (x-G)p 2 . (22b) 
Accordingly, both Hq and Hq 2 '* contribute to the pole part of I: 

w+o-n/o _ u _ A r x /•1-x+G 

1 = " T^n^nS L d ^ L dG L d y( H !) + Hi 2) )" +ff " 4 + Finite Terms , (23) 

or, since 



( 2vr )2(c+a) J Q J Q J Q 



Hence, 



dy(...)] ^ = [/ dy(...)] (24) 



I = 2F(4 u a) f 1 d xx- CT / " dGG w+CT - 3 + F.T. , (20) 



(4tt) 2 Jo " Jo 
2f(4 - u - a) + ^ _ /•! dx ^_ 2 + p T ^ (27) 



o 



or, finally 



div 



(47r) 2 p 2 p 2 

d 4 q 1 -4r(2 -uo-a 



(27r) 4 q 2 (q + p) 2 p 2 (q + p) 2 p 2 p 2 (4vr) 



<7-»(l/2) + 



U)^(3/2)H 



(28) 



which is the result previously given in Table A.l of Ref. Notice that the above noncovariant 
integral has both a three-dimensional and a four-dimensional nonlocality. 

4.2 The five-propagator integral 

We now turn our attention to the massive five-propagator integral I in Euclidean space, 

I=/J*4_ ® __ ,k = p'-p, i = 1,2,3 , (29) 

J (27r) 4 q 2( q _ k )2[( p + q )2 + m 2]^2(^ _ k )2 

which is seen to possess two noncovariant and three covariant propagators. Employing a special 
parametrization we may re- write I as 



{fl n rl— x+G /-x— G ~) noo f rl 4 n 

/o dX i dG /c dy i dV j/o dAA /(2^" e < 3 °> 



where: 



E = Ax(q 4 + r 4 ) 2 + A(q + r) 2 + AZ , 



vpi-yk;, r 4 = x 1 (vp 4 - Gk 4 ) , i = 1, 2, 3 , (31) 



Z = v[p 2 + m 2 + 2Gp 4 k 4 /x - vp 2 + 2yp • k] + y(l - y)k 2 + G(x - G)k 2 /x - v 2 p 2 /x , 
so that 



I = {...} jH dAAV AZ J d 3 q S e~ A( ^ 2 J dq 4 e 



Ax(q 4 +r 4 ) 2 



Defining 

Qi = $ + Q 4 = g 4 + r 4 , i = 1, 2, 3 , 

and implementing split dimensional regularization, we find that 



I={ ... } /; d AAV-/^( Ql - ri )e-^/ M 



4 e -AxQ| 



2ct 



oo 

4 — uj — tj ( _ \, — o- — AZ 



= {-} (2^)/„ dAA (— ri)x -<7 e , (32) 

—Yi^ — (i! — rr\ r 1 f x /-1-x+G /-x-G 

= i d -"'/o dG /o *i dv(vp,- ykl )Z— , (33) 

where we have made use of the following formulas: 

J d 2w Q exp(-AQ 2 ) = t^A"" , A > ; 

J d 2w Q Qi exp(-AQ 2 ) = , i = 1, 2, 3 ; 

J d 2CT Q 4 exp(-AxQj) = tt ct (Ax)- ct ; (34) 
/ dAA 4 """" exp(-AZ) = T(5 - u - <t)Z w+ct - 5 . 

We re-iterate that Q is defined over 2u;-dimensional complex space, and Q 4 over 2<r-dimensional 
complex space. 

Our next task is to extract the pole part from the integral in Eq. (33), by locating the 
singularities of the integrand in four-dimensional parameter space. Proceeding as in Section 
4.1, we see that Z in Eq. (31) vanishes for the following sets of integration parameters: 

Case 1: 

Z — > 0, if v = 0, y = 0, G = 0, in which case 

Z -> Z^ = v(p 2 + m 2 ) + yk 2 + Gk 2 . (35) 
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Case 2: 

Z — > 0, if v = 0, y = 1, G = x, in which case 

Z -> z{ 2) =v(p 2 + m 2 + 2p-k) + (l-y)k 2 + (x-G)k 2 . (36) 

To obtain Zq 2 -' in Eq. (36), it is convenient to introduce new parameters Y = 1 — y, and 
t = x — G, so that 

Z^O, if v = 0,Y = 0,t = , (37a) 

and 

if = v(p 2 + m 2 ) + Yk 2 + tk 2 , p 2 = p 2 + 2p • k , (37b) 
which is analogous to in Eq. (35). 
Accordingly, the total expression for I in Eq. (33) reads: 

_r(5 - u - a) r 1 r /-i-x+g /-x-g 

.[( Z ( 1 ))^-5 + ( Z (2))o;+a-5] _ (3g) 

The divergent part of the contribution from Z^, at v = 0, y = 0, G = 0, can be shown to 
vanish. 

To derive the pole contribution from Zq 2 - 1 , we replace the parameters {x,G,y,v} by {x,t,Y,v}, 
whence 

pi pX pi — x+G px—G pi px pt pt 

/ dx / dG / dy / dv = - / dx / dt / dY / dv , (39) 

Jo Jo Jg Jo Jo Jo Ji+t-x Jo 

and 

I = r ^~^~ (7) ^ dxX-^ f dt r dY t dv( VPi + Yk; - ki) 

(4tt) w + ct Jo Jo Ji+t-x Jo 

• [v(p 2 + m 2 ) + Yk 2 + tk 2 ]^" 5 . (40) 

Since the integration proportional to (vp; + Yk ; ) leads to a finite value, Eq. (40) reduces to 
the expression 



(4tt) w + ct ) v ; Jo Jo Jl+t-x 

dv 

[v+(p 2 + m 2 )- 1 (Yk 2 + tk 2 )] 



r x dv 

/ , (41) 

Jo f v + fr> 2 + T n 2N - i rYk 2 + tk 2 n 5 -"- CT 



or, eventually, to 



div 



/ 



d 4 q ft 



(27r) 4 q 2 (q - k) 2 [(g + p) 2 + m 2 ]q2(q - k) 2 
— 2kil* 



k 2 k 2 [(p + k) 2 + m 2 ] 



,<7->(l/2)+ W ->(3/2)+, (42) 



where I* = T(2 — uj — a)/ (47r) 2 . The above integral is seen to be nonlocal in both k 2 and k 2 . 

The method described here is applicable to all four- and five-propagator cases. Particularly 
challenging is the basic five-propagator integral, whose pole part consists of two distinct terms, 
namely 



div 



d 4 q 

) 2 + m 2 ]q 2 (q — k) 2 

(43) 



(2vr) 4 q 2 (q - k) 2 [(p + k) 2 + m 2 ]q 2 (q- k) 2 
-21* / 1 1 \ 



k 2 k 2 \P 2 + m2 (p + k) 2 + m 2 y 

5 BRST identity in the Coulomb gauge 

It remains to convince ourselves that S(p) and A M in Eqs. (10) and (16) satisfy the appropriate 
one-loop BRST identity. In order to do that, however, we have to know what the 'appropriate 
identity' really is. The expression usually quoted in the literature has the structure 

(p' - pf A>', p) = -S(p') + S(p) , (44) 

and the question is: does this identity also hold in the Coulomb gauge, where ghosts are known 
to play a significant role? The answer to this question, surprisingly, is no. 

According to Taylor, the correct BRST identity also involves ghost contributions, Gi(p',p) 
and G2(p',p), and is given by |15|, [16 



(p' - P )^(p',p) + Gi(p',p) + G 2 (p',p) = -S(p') + E(p) . (45) 

The corresponding ghost diagrams are depicted in Figs. 4a and 4b, with solid lines denoting 
quarks of mass m, curly lines gauge bosons, and with broken lines representing ghosts; q is the 
internal momentum, as before, and p,p' are external momenta. It remains to compute the 
functions Gi, G 2 . 

The ghost diagram depicted in Fig. 4a leads in Minkowski space [[TJJ to the expression 
Gi( P ',p), 



Gi(p',p) = [factors] J 



d 4 q i L 5 .. + (q- k )i(q- k )j 



(2tt) 4 [(q - k) 2 + is] 



(q-k) 2 



•^o-gi7j //, , — ^(^-m),£>0,k = p' -p . (46) 

Applying the technique of split dimensional regularization to all Coulomb-gauge integrals, and 
utilizing the formulas in the Appendix, we find that the divergent part of Gi is zero, i.e. 

divGi(p,p) = . (47a) 
A similar result may be established from Fig. 4b for G 2 : 

divG 2 (p / ,p) = . (47b) 
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In summary, neither Gi nor G2 contributes to the BRST identity in Eq. (45), which is 
obviously satisfied by S(p) and A^(p',p). After completion of this work, we discovered that an 
identity closely resembling the identity in Eq. (45) had been studied by Muzinich and Paige 
|I~T|1 and, more recently, by Newton II7J using a new approach. 



6 Genuine energy integrals at two loops 

The purpose of this section is to discuss briefly the application of split dimensional regularization 
at the two- loop level |Tj| [19|]. In particular, we would like to know to what extent the new 
Coulomb-gauge technique is capable of regulating both the divergences from the notorious 
two-loop energy integrals || || |7j, ||, as well as the ordinary UV divergences. The safest 
way of getting at least a partial answer to this question is to embark on an explicit two-loop 
calculation and then see what happens to the ambiguous energy integrals. A program of this 
kind has recently been initiated by the author, but the final results won't be available for at 
least a year or two. In the meantime, we shall confine ourselves to a few general remarks. 

Our first comment deals with the measure d 4 q d 4 k of a two-loop integral. Applying rule (1) 
for one-dimensional integrals, we employ exactly two complex-dimensional parameters uj and 
a, thus replacing the measure by 

d 2(^+ CT ) q d 2(-+a) k = d 2 CTq() d 2^ d 2<x ko d 2^£ . ( 4 g) 

the limits a — > (l/2) + , and uj — ► (3/2) + are taken after all integrations have been executed. 

Next, a general two-loop Coulomb-gauge integral is expected to give rise to simple and 
double poles, proportional to T(2 — (uj + a)) and [r(2 — (u> + o - ))] 2 , respectively; notice that uj 
and a appear additively. 

We now turn our attention to a typical energy integral, such as [0] 

.dqo, dko ko qo (4Q) 

J 2tt J 2tt (kg _ k 2 + ie) (qjj - q 2 + ie) 1 ' 

which occurs for the first time at two loops. According to Doust and Taylor |], |^, || , standard 
dimensional regularization is incapable of handling the ambiguities in integrals like (49). In 
order to treat (49), Doust and Taylor use a special regularization, called ^-regularization, which 
assigns to (49) a particular value, such as 1/12. Our approach to this problem is somewhat 
different: by applying split dimensional regularization to a certain two-loop integral, we hope to 
regulate both its ordinary UV divergences, as well as the ambiguities from its energy integrals. 
However, before discussing a specific example, we shall briefly indicate how energy integrals 
such as (49) arise in practice. 

Let us consider the two-loop Yang-Mills self-energy in the Coulomb gauge, specifically the 
sunset diagram (it possesses overlapping divergences). The amplitude for this challenging dia- 
gram reads as follows (we employ the notation of Ref. [2"0|): 



n£,(p) = - // Vf4 GfM G%( V - k - q) Gg(k) Vg* , (50) 
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where G^(q) represents the gluon propagator in the Coulomb gauge, Eq. (4), and 



v?4 = -ig 2 [f abc f cdc (v^p-^ P ^) 



+ P ce i hde (5 pX 5 ap - S w S Xa ) + f adc f cbc ( 5,J Xp - 8^5 ap )\ 



adcrcbc/ 



(51) 



denotes the four-gluon vertex of zero-loop order. Reduction of the integrand yields, in Euclidean 

space, 



where: 



n£(p) = (coeff.HV//. 



d 4 q d 4 k 



-M 



q 2 k 2 (p — k — q) 2 
x [A paA5 (q,k;n) + B paA5 (q,k;n)]C (T/3 (q,k;p,n) , 



(52) 



+ 35 pX 5 ap (25 a /?5^ — 5 a (_5(3u — Savbpe) ; 



(53) 



A pa Ae(q,k;n) = 5 pa 5 Xi . - 



k • n S pa n x k^ k • n 5 pa n^k x q • n 5 x/ :n p q a 



k 2 



k 2 



q 2 



q-n^n Q q p q • nk • n (n p q a + n a q p )(n A k f + n ? k ; 



q 2 



+ 



cpk 2 



C CT/ 3 = 5 a /3 + 



'^R^R^ - R • n (R CT ng + R^n CT )' 
, ^R 2 " , 



R<7 = (p k q) tj 



(54) 



(55) 



Bp a Ag consists of four terms, but is not required for our purposes. Substitution of expressions 
(53)- (55) into Eq. (52) leads to three kinds of integrals: (i) covariant integrals, (ii) noncovariant 
Coulomb-gauge integrals, and finally, (iii) noncovariant Coulomb-gauge energy integrals. At 
present, we are only interested in extracting a few integrals in category (iii). 

To locate a ko energy integral, for instance, we multiply the seventh term in M paXpa ^j3 U by 
the second term in A paX ^, and then contract the resulting expression with the n 2 R (7 R / g-term in 
C a/3 . Thus, 



TpJ = (65pA5<rp(W%/) 



-k ■ n S pa n x k^ 
k 2 



n 2 R a R 



-R 2 



6 k • n R 2 n M k^ 

" M 2 ' 

so that the corresponding integral from Eq. (52) reads 



(56) 



I(i) 



n, 



d 4 q d 4 k k • n k u 
q 2 k 2 k 2 (p-k-q) 2 



(57) 
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or, in Minkowski-space, 



t(D ( \ f d4q f d ^ f dk ° k ° k » 1 1 c« 



leading to the ko energy integral 



dk k 



(59) 



eg - k 2 + ie) 

Similarly, multiplying the fourth term in M^xpo^pv by the fourth term in A pa \^, and contracting 
lp in C a /3, 



with n 2 R a RB in CLg, we find that 



T$ = 6q-nRVW(q 2 F?) • (60) 
The corresponding integral contains a go energy integral of the form 

dq qo 



(Oo - q 2 + ie ) 



(61) 



To conclude this section, we shall illustrate how split dimensional regularization handles 
energy integrals like (49), 

rd^dko ko qo 

7 2tt J 2vr (kg _ k 2 + i e ) (qjj - q 2 + ie) K J 

which corresponds to the integral (1.3) in Doust [/j. Defining first fco and go each over 2a- 
dimensional space (and labeling them and Q^, respectively), and then defining k and q 
each over 2uj- dimensional space (labeling them K and Q, respectively), we may rewrite (62) as 



d 2a Q r d 2a K K u Q u 

-^^r / 7 — \o £ , e > 0; 63 

27T 2 - 7 27T 2 - (#2 _ K 2 + i£ ) (Q2 _ Q2 + 



if 2 and Q 2 denote the squares of the 2<r-dimensional vectors and Q^, respectively, while K 2 
and Q 2 are the squares of the 2a;-dimensional vectors K and Q, i.e. the squares of the original 
3- vectors k and q. Wick-rotating l M to Euclidean space, we obtain 

= ^^Q_ /• rf 2 -K ^ Q M 
E J (2vr) 2 - 7 (2tt) 2 - (K 2 + K 2 ) (Q 2 + Q 2 ) ' 1 J 

Since each one of the integrals in (64) is now well defined, we may apply standard dimensional 
regularization to get 

d 2a Q Q„ n f d 2 °K K„ 



(27r) 2 -(Q2 + g 2) °' J (2n) 2 ° (K 2 + K 2 ) °' (65) 

so that 1 = 0. Consequently, in split dimensional regularization, the ambiguous energy integral 
(62), or (49), has the value zero. 
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7 Conclusion 



In this paper we have carried out a second test of split dimensional regularization by evaluating 
to one loop the pole parts of the quark self-energy E(p), and the quark-quark-gluon vertex 
A M (p', p). The technique of split dimensional regularization was specifically designed to regulate 
Feynman integrals in the noncovariant Coulomb gauge y ■ A a (x) = 0. Its most important 
single feature is the use of two complex parameters, u and a, which permit us to regulate the 
divergences of certain energy-integrals. 

The results in this paper may be summarized as follows. 

1. Split dimensional regularization enables us to compute unambiguously all relevant one- 
loop integrals in the Coulomb gauge. This comment applies especially to the noncovariant 
four- and five-propagator integrals appearing in the vertex function A At (p / , p) (Section 4). 

2. The conventional BRST identity in Eq. (44) should by rights be replaced by the general 
BRST identity in Eq. (45). The latter contains two new functions, Gi and G2, which 
correspond to the ghost vertex diagrams shown in Figs. 4a and 4b, respectively. An 
explicit calculation of Gi(p',p) and G 2 (p',p) reveals, however, that the pole parts of 
these functions are actually zero. 

3. Computation of T^(p',p) and r^(p',p) involves several complicated Coulomb-gauge inte- 
grals, some of which turn out to be nonlocal. Nevertheless, both and are strictly 
local and respect, together with S(p), the correct BRST identity in Eq. (45). 

4. Although ghosts generally play an important role in the Coulomb gauge, they do not 
contribute explicitly to the BRST identity relating A M and S(p). 
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Appendix 



The following Coulomb-gauge Feynman integrals, defined over 2(u + <r)-dimensional complex 
Euclidean space, are useful in the evaluation of S(p), and r 2 in Eqs. (10), (13) and (15), 
respectively. We employ the abbreviations 



and 



dq = d d qdq 4 /(27r) 4 , 



T(2-w-a) 
(4^)2 



<r-(l/2)+ 
w-»(3/2)+ 



div 
div 
div 
div 
div 



dq g^q 4 



q 2 [(q + p) 2 + m^q 2 6 



dq ql 



q 2 [(q + p) 2 + m 2 ]q 2 2 



/fp' 



dq 



(p' + q) 2 [(q + p) 2 + m 2 ]q 2 
dq ftq 4 



div 



(p' + q) 2 [(p + q) 2 + m 2 ]q 2 
dq 



, 

0,i = l,2,3, 



div 
div 
div 
div 
div 
div 



[(p' + q) 2 + ni 2 ][(p + q) 2 + m^q 2 
dq 



J \(V> 



0, 
1 



[(p' + q) 2 + m 2 ] [(p + q) 2 + m 2 ]^ 3 



dq 



-21* 



q 2 (p' + q) 2 [(p + q) 2 + m^q 2 (p') 2 (p 2 + m2 ) 

d q gigj = g 

q 2 (p' + q) 2 [(p + q) 2 + m^q 2 

d q gigj = q 

q 2 [(p' + q) 2 + m 2 ][(p + q) 2 + m^q 2 



dq 



-21* 



/ 



q 2 [(p + q) 2 + m 2 ](q + p'^q 2 (p') 2 (p 2 + m 2 ) 
dq q { 



q 2 [(p + q) 2 + m 2 ](q + p') 2 q 2 



0, 
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diy f dq gjq 4 = 

J q 2 [(p + q) 2 + m 2 ](q + p') 2 cp 

div r d q m _ -spfcj' 1 * 



q 2 (p' + q) 2 [(p + q) 2 + m 2 ](q + p') 2 ^ 2 (p') 2 (p') 2 [(p - p') 2 + m 2 ] ' 
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Figure 1: Quark self-energy graph. Solid lines denote quarks, curly lines gluons. 




Figure 2: QED-like quark-quark-gluon vertex graph. 




Figure 3: Non-Abelian quark-quark-gluon vertex graph. 
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Figure 4: (a) Ghost contribution ||T5| to the IPI three-point function, with broken lines repre- 
senting ghost particles, (b) Mirror image of the ghost diagram depicted in Fig. 4a. 
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